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ABSTRACT 

We present a survey of the X-ray emitting ejecta in the Cassiopeia A supernova remnant based 
on an extensive analysis of over 6000 spectral regions extracted on 2.5-10" angular scales using the 
Chandra 1 Ms observation. We interpret these results in the context of hydrodynamical models for 
the evolution of the remnant. The distributions of fitted temperature and ionization age, and the 
implied mass coordinates, are highly peaked and suggest that the ejecta were subjected to multiple 
secondary shocks following reverse shock interaction with ejecta inhomogeneities. Based on the fitted 
emission measure and element abundances, and an estimate of the emitting volume, we derive masses 
for the X-ray emitting ejecta and also show the distribution of the mass of various elements over the 
remnant. An upper limit to the total shocked Fe mass visible in X-rays appears to be roughly 0.13 
M0, which accounts for nearly all of the mass expected in Fe ejecta. We find two populations of Fe 
ejecta, that associated with normal Si-burning and that possibly associated with a-rich freeze-out, 
with a mass ratio of approximately 2:1. Essentially all of the observed Fe (both components) lies well 
outside the central regions of the SNR, possibly having been ejected by hydrodynamic instabilities 
during the explosion. We discuss this, and its implications for the neutron star kick. 

Subject headings: hydrodynamics - ISM: individual (Cassiopeia A) - supernova remnants - X-rays: 
ISM 



1. INTRODUCTION 

Theory and observations are merging toward consensus 
that core-collapse supernova explosions are intrinsically 
asymmetric, even if the progenitor had an initially sym- 
metric configuration. Recent three-dimensional simula- 
tions show that strong dynamical interactions between 
the various burning shells can lead to large asymme- 
tries in the proge nitors, independent of the symmetry 
of the explosion ([Arnett fe Meakinl [20lTI ) . Nor is the 
explosion itself li kely to be spherical, with conve ctive 
instabilities (e.g. iHerant et al| [l992; Herant 199 5]) and 
insta b ilities at the accreti on shock (e.g. iBlondin et all 
120031: iFoglizzo et all 120071 ) all possibly acting to cre- 
ate asymmetries. It is further anticipated that such 
asymmetries may hold t he key to the co r e-collapse ex- 
plosio n mechanism (e.g. lOtt et all [2008; Br andt et al.l 
1201 If ). They will be manifested not only in the dis- 
tributions of the ejecta mass and velocity, but also in 
the recoil of the nascent neutron star, which is ob- 
served wit h measured space velocit i es of up to 1000- 
1500 km/s (lArzoumanian et al ll200l iHobbs et al.ll2005l : 
iFaucher-Giguere fe Kas pi 2006). Simulations show that 
a neutron star kick that has a purely hydrodynamic ori- 
gin can balance the total momentum of the anisotropic 
component, and is expecte d in the direction oppos i te the 
fastest shock exp ansion ([Wongwathanarat et al.l 120101 : 
IScheck et al][200l . 

Observations of core-collapse supernovae indicate a sig- 
nificant degree of explosion asymmetry via measurements 
such as ti me-dependent polarization changes or Do ppler 
velocities ([Mazzali et alll2006l; iLeonard et al.ll2006D . but 
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explosion asymmetries are also revealed in sufficiently 
young supernova remnants (SNRs) through their shock- 
heated ejecta. The most detailed view so far available of 
the ejecta in a young core-collapse SNR is provided by 
Cassiopeia A, which at 330 years is the penultimate of the 
known Galactic supernova remnants. Its reverse shock 
has already progressed deeply into the explosively pro- 
duced nucleosynthesis products, aided by the strong pre- 
supernova mass loss incurred by its proge nitor through 
the l ikely action of a binary companion ([Young et al.l 
I2006T ). 

That Cas A was produced by an asymmetric explo- 
sion is by now well-established, and most recently by 
light echo observations which sample the explosion hun- 
dreds of years after the fact in different directions. These 
show variations in the ejecta velocities of ~ 4000 km/s 
(Rest et al. 2010). Chandra X-ray observations have 
also played an important role in revealing the complex- 
ity of Cas A's X-ray emitting ejecta. The X-ray emitting 
Si ejecta show a bipolar structure with jet-like features 
(|Hwang et al.ll2004l:lvink et al.ll2004|:lLaming et al.ll2006l) 
similar to that seen in optical (jFescn 2001, a nd refe rences 
therein) and infrared emission (jHines et al.l[200l . The 
X-ray ejecta spectra reveal the imprint of roughly a factor 
of two asymmetries i n the deposition of explo sion energy 
around the remnant ([Laming fe Hwang 112003! ) . similar to 
the distribution of kin etic energies and position s of fast 
moving optical knots (Ha mmell fe Fe sen 2008). Early 
Chandra observations identified regions dominated by 
emission from Fe ejecta well outside the projected loca- 
tion of Si emission that were proposed to b e sites of over- 
turn between the nucleosynthesis layers ([Hughes et al.1 
2000). This conclusion is apparently supported by the 
X-ray Doppler measu rements showing highe r velocities 
of Fe compared to Si ([Willingale et all 120021 ). More re- 
cently, the dynamical structure of the ejecta has been 
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inferred in remarkable detail by the multi-wavelength 
(X -ray, optical, in f rared ) Doppler-shift maps compiled 
by iDeLanev et all (|2010T ) . These show the inner ejecta 
to be unshocked and in a flattened distribution ([Si II]), 
outlying matter to be arranged in rings on a spherical 
surface ([Ar II], [Ne II]) , and outflows emerging out of 
the surface of the remnant (Fe ejecta, jet, and counter- 
jet). Their detailed analysis shows that the outflows of 
X-ray emitting Fe ejecta are not material mixed outwards 
from overturning of the ejecta layers, but rather, mate- 
rial that has pushed through the overlying ejecta and is 
now encircled by rings of material from outlying nucle- 
osynthesis layers. 

Large-scale studies of the X- ray emission inCas A 
have already be e n und ertaken bv lStage et al.l ()2006l ) and 
IHelder fe Vink I (]2008[ ) , but these focus on the properties 
of the nonthermal emission. The most comprehensive ac- 
counting of the X-r ay emitting ejecta in Ca s A to date 
was carried out by fWillingale et al.l (|2003l ). who used 
XMM-Newton CCD observation s to examine 225 region s 
of fixed 20" size across Cas A (jWillingale et alj|2002ft . 
They fitted a two-component thermal model through- 
out to represent hot shocked circumstellar material plus 
cooler shocked ejecta, and inferred a total ejecta mass of 
2.2 M Q . They interpret the Fe K emission as forming a 
bipolar double cone that is associated with ejecta bullets 
that have broken beyond the forward shock. 

A detailed accounting of the X-ray emitting ejecta on 
fine angular scales was one of the primary goals of a 1 Ms 
Chandra observation of Cassiopeia A that was obtained 
in 2004, and is the subject of this paper. Our spectral 
survey includes over 6000 spectral regions extracted on 
2.5-10" angular scales. It is distinguished from previous 
Chandra studies of the ejecta in Cas A, which have relied 
either on spectral imaging or the detailed spectral anal- 
ysis of a relatively small number of regions. We are able 
to cover a larger area of the remnant with much higher 
sensitivity and angular resolution than was available to 
iWillingale et al.l ((2002). In particular, we are able to 
show the distribution of the X-ray emitting Fe ejecta in 
detail. The Fe ejecta, being produced just outside the 
collapsing core of the supernova, provide valuable clues 
to understand important details of the explosion, includ- 
ing the recoil of the ejecta with the neutron star. 

We interpret the spectra in the context of the hy- 
dro dynamical models applied by us previously to Cas 
A (lLaming fe Hwang I 120031: IHwang fe Laming"! 120031 : 
lLaming et al.l 120061 : iHwang fe Laming 1 120091 ). Cas A 
particularly lends itself to qua ntitative analysis by virtue 
of its well-const rained age dThorstensen et all 120011 : 
iFesen et al.l 12006ft and distance (|Reed et all 119951) . Re- 
cent one-dimensional hydrodynamical models for the 
interaction of the remnant with a circumstellar wind 
medium indicate that the reverse shock in Cas A has al- 
ready interacted wi th a si gnificant fraction of the ej ecta 
(jChevalier fe Oishi|[200l . lLaming fe Hwang I (|200l ap- 
ply their models directly to Chandra X-ray spectra, and 
also infer that there is very little of the unshocked ejecta 
remaining. Thus a large fraction of the ejecta in Cas 
A are accessible to X-ray observations. The exceptions 
are the small amount of ejecta that have not yet been 
shocked, ejecta at very low densities (which we will dis- 
cuss later), any ejecta that might have been subject to 
rapid cooling by thermal instability, and and ejecta of 



very high density that have been shocked to lower tem- 
peratures and thus em it at longer wavelengths as optical 
knots (e.g.|Ffesen||200lD. 

The organization of the paper is as follows: Section 2 
gives an overview of the analysis procedure, with expla- 
nations as to the spectral grid and background, and most 
importantly how the spectra were organized according to 
the fitting model required. We then proceed to discuss 
the results, from maps for the remnant as a whole to a 
few observations about properties of the forward shock- 
dominated regions, and then of the ejecta. Section 3 dis- 
cusses issues concerning the models and the refinements 
that have been made for this work compared to previous 
work. Section 4 then presents the results of the mass cal- 
culation and issues connected to it, such as the presence 
of unshocked ejecta and implications for the neutron star 
kick. 

2. ANALYSIS PROCEDURE 

We use the 1 Ms Advanced CCD Imaging Spectrom- 
eter (ACIS) observation obtained by the Chandra X- 
ray Observatory in February and May 2004 (one of the 
first Chandra VLP, or Very Lon g Project, observati ons). 
This observation is described bv IHwang et al.1 (|2004[ ). Its 
most salient features are that it was obtained in nine ob- 
servation segments (OBSIDS) and in GRADED mode, 
wherein the CCD detection events are characterized on- 
board the spacecraft before telemetry. It was necessary 
to use GRADED mode to reduce the telemetry load, 
due to the high source count rate (the observation ac- 
cumulated some 280 million photon events in 980 ks). 
The main disadvantage of using GRADED mode is that 
detailed corrections of detector problems such as pulse 
pile up are not possible. The data were processed with 
CALDB Version 3.2.2. 

In this section, we outline the analysis procedure, 
which includes definition of the spectral grid, consider- 
ations for the subtraction of the background, the classi- 
fication of forward shock and ejecta-dominated spectra, 
the selection of appropriate spectral models, and the re- 
sults of the fitting. 

2.1. Spectral Grid 

To carry out our spectral survey of the ejecta in Cas 
A, we define a grid of 10" boxes across the entire super- 
nova remnant (SNR), optionally subdividing by factors 
of two down to 2.5" depending on the number of counts. 
The regions are thus square or rectangular with sides of 
length 2.5, 5, or 10". We excluded the regions dominated 
by the central compact object, and a few other regions 
with low numbers of counts. This gave us 6202 spectra 
for analysis. Although this is a very large number, the 
data would have readily allowed at least a factor of two 
finer grid sizes down to nearly the Chandra angular res- 
olution in some regions. This being the first attempt at 
an analysis of the complex thermal emission on such a 
fine angular scale, we tried to keep the number of grid 
regions relatively tractable. 

Each spectrum was then extracted individually from 
all nine OBSIDS comprising the Ms observation, and the 
corresponding individual photon redistribution matrices 
and effective area files computed. The individual spec- 
tra were then added, and the individual response files 
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weighted and added, to obtain the final spectral and re- 
sponse files for each grid region. The pulse-height spectra 
were further binned to yield a minimum of 25 counts per 
bin to apply Gaussian statistics. 

The filamentary and knotty characteristics of Cas A's 
X-ray emission, its high surface brightness, and the ex- 
cellence of the Chandra mirrors, result in pulse pile-up 
at the brightest ejecta knots. Pulse pile-up can pro- 
duce spurious line features, for example, with energies 
that correspond roughly to the summed energy of the Si 
and/or S He a blends that are so prominent in Cas A. 
Pile-up of continuum photons will also skew the spectral 
shape and result in higher fitted temperatures. 

A detailed correction for pulse pile-up cannot be per- 
formed for these GRADED mode data, and we made no 
attempt to account for these effects for the thousands 
of spectra examined here. We can, however, estimate 
qualitatively the number of our spectra for which pile-up 
is a significant problem. For the vpshock fits described 
later (in section 2.3), we consider the distribution of the 
fitted temperature against the counts per pixel for each 
spectrum. The average fitted temperature shows a clear 
trend of increasing for spectra with more than about 6000 
counts/pixelin980kfl For the 25 regions with the high- 
est count /pixel ratios (> 6000 counts/pixel), the average 
fitted temperature is 2.25 keV, compared to an average 
of 1.7 keV for the remaining spectra. Inspection of these 
spectra show the problems in fitting the line emission de- 
scribed above; these problems are not apparent for lower 
count-rate spectra. Since the pile-up occurs noticeably in 
only a limited number of regions (~25), we expect that 
its uncorrected effect on the total mass is not large. 

2.2. Background 

The off-source regions outside the SNR show a line- 
rich spectrum similar to that of the SNR itself. The 
background is thus dominated by the bright source spec- 
trum due to scattering by intervening dust, the CCD 
readout, and to a lesser extent, the wings of the mir- 
ror point spread function. Dust scattering has been 
shown to be significant f or Cas A based on RQSAT and 
Einstein observations bv IPredehl fc Schmittl ()1995f ) and 
iMauche fe Gorensteinl (|1989D . 

Since the source spectrum, and hence the background, 
varies both in form and brightness around the remnant, 
we extract 16 separate background spectra by subdivid- 
ing a 0.4' thick, ~ 3.5' radius circular shell surrounding 
the remnant so that each segment covers roughly 20 de- 
grees in azimuth. We then assign a background spectrum 
for each source spectrum based on the source region's 
azimuthal position angle. Where the source extends be- 
yond the background extraction radius, near the north- 
east "jet" and its southwest counterpart, we use the sum 
of the two nearest background regions to either side. 

We choose to subtract the background from the source 
spectrum rather than to fit a model to each background 
spectrum. A model for the background would have to 
include a thermal component for the scattered source 

3 The count /pixel cut-off we discuss here is of course specific 
for this particular observation. For reference, the brightest pixels 
contain on the order of 16000 counts, corresponding to a count 
rate of 0.05 counts per 3.24 s frame. This corresponds roughly to 
pileup fractions of less than a few percent according to the Chandra 
Proposers' Observatory Guide. 



spectrum, whose normalization would inevitably be dif- 
ficult to constrain. The normalization for the particle 
background component is well-constrained by the spec- 
trum at high energies, but the effect of scaling the scat- 
tered thermal background in the same way as the parti- 
cle background is similar to simply subtracting the total 
background. 

Simple subtraction of the off-source background means 
that the background is generally underestimated, and by 
different amounts depending on the photon energy and 
the dominant scattering process. The precise determina- 
tion of the background is thus a sophisticated problem 
that is beyond the scope of our study. The ejecta survey 
which is the focus of our efforts will be less sensitive to 
the precise background subtraction because most of the 
ejecta regions are bright, with strong line emission; the 
subtleties of the background subtraction are less impor- 
tant here than for inferences about the origin of weak 
lines. 

2.3. Spectral Classification 
2.3.1. Spectral Map 

To provide a consistent spectral characterization of 
the entire remnant, we fit to every spectrum in our 
grid a single-component plane-parallel shock model with 
variable element abundances (XSPEC model vpshock; 
Borkowski et al. 2001) modified by interstellar ab- 
sorption (XSPEC model wabs; Morrison & McCammon, 
1983). The vpshock model incorporates ionization ages 
from zero up to a fitted maximum value. In previou s 
work (|Laming fe Hwang 12003 iHwang fe Laming 12 003). 
we had used simple nonequilibrium ionization (NEI) 
models characterized by a single temperature and ion- 
ization age, but there the extraction regions were gener- 
ally about 3" and the exposure time only 50 ks. We 
find that for the much longer exposure time of these 
data, and the sometimes larger extraction regions, such 
models are no longer acceptable. Because most of the 
regions have ejecta-dominated spectra, we have not in- 
cluded any elements lighter than O in the model, tak- 
ing O to b e the prima r y sou r ce of the continuum in the 
man ner of iVink et al.l dl996H. iLaming fe Hwang"! (j2003[ ) 
and iHwang fc Laming I ()2003l ) . The abundances of ele- 
ments heavier than O are generally varied (with Ni linked 
to Fe, Ca to Ar, and sometimes Ar to S), while O is 
he ld fixed at the solar value, using the solar abundances 
of lAnders fc Grevessd (|1989l fl. The redshift is allowed 
to vary to either positive or negative values, given that 
the ejecta in Cas A have significant bulk motions (e.g., 
iMarkert et aTHl983tlWillirigale et alj|2002ft . and is gener- 
ally driven by the Fe L emission and the strong Si blend. 
Given that there may be calibration problems near Si 
(see the discussion in DeLaney et al. 2010), we consider 
the fitting of the redshift more as an aid to achieving 
good fits rather than a reliable measure of the actual 
line-of-sight velocity. The fitted redshift does, however, 
reproduce the gross features of the bulk motion as mea- 
sured by other means. A Gaussian smoothing scale for 

4 In practice, the emission is calculated as if all elements, specif- 
ically H, were included at their solar abundance ratios. The indi- 
vidual contribution of each element is then further weighted by the 
abundance value relative to solar, so that those elements with zero 
abundance values do not contribute. 



4 



Hwang & Laming 



the spectrum was also fitted, given that the observed 
strong spectral lines are generally broad because of bulk 
motions. 

While we do not present errors for the fitted parame- 
ters here, these were computed for key parameters. We 
have found that this is a necessary step to optimize find- 
ing the true minima: the x 2 terrain is generally rather 
rugged for these fits. The redshift and ionization age 
in particular are prone to settle at secondary minima, 
as might be expected, considering that these parameters 
both have a strong effect on the energies of the strong 
line features. 

Here we note some of the general features of these re- 
sults based on images shown in Figure f, where fitted 
model parameters such as temperature, ionization age, 
and various element abundances are shown for each spec- 
tral region. Histograms and a scatter plot for the fitted 
temperatures and ionization ages are shown in Figure 
2. The column density increases systematically to the 
west by about a factor of 2 or more and reaches a max- 
imum at the western extremity of the remnant. This 
is entirely consistent with pr evious work ( Kcohan e et al.l 
19961 iWillingale et all 12002 ). There are also smaller re- 
gions in the center with significant localized column den- 
sity enhancements of about 50%. The temperature dis- 
tribution has two peaks, which are seen most clearly 
in the histogram in Figure 2. The temperature map 
shows that the lower peak is associated with ejecta- 
dominated, line-emitting regions and the higher with 
the forward shock-do minated regions that were identi- 
fied in previous work dGotthelf et alj [200lt iStage et al.l 
[20061 IHelder fe VinYl 120081 ) . (We will discuss the for- 
ward shock in the following section.) The ionization age, 
by contrast, is remarkably narrowly distributed overall, 
though coherent regions with high ionization ages can be 
identified. In the western half of the remnant, these high 
ionization regions appear to be associated with the for- 
ward shock. The highest ionization ages are found in the 
east, however, and there they are clearly associated with 
Fe-enriched ejecta, as is apparent from comparing the 
ionization age and Fe abundance maps shown in Figure 
1. 

We will focus our attention primarily on the abun- 
dance maps for Fe and Si, although such maps are also 
shown for all the elements for which we fitted the abun- 
dances. The maps for S and Ar are on the whole sim- 
ilar to that of Si. The abundance maps for Si and 
Fe are distinctly different, but both feature three main 
lobes of ejecta located north, east, and west. The Si 
is also extended to the northeast along the northeast 
ejecta "jet" , while Fe is particularly distinct in the east- 
ern region compared to the Si, and shows its highest 
enhancement in the outermost parts of this region. In 
general, the abundance maps we obtain are strongly 
reminiscent of the cor r espon ding l ine images shown b y 
iHwang. Holt, fe Petrel pOOOl ) and IHwang et all (|200l . 
and also r esemble the less detaile d abundance maps ob- 
tained by IWHIiniakitaD J2QQ3) • Those authors had 
previously noted that the line emission and element 
abundance patterns for Si, S, Ar, and Ca in Cas A are 
similar to each other. 

The Ne and Mg maps have a distinct character in that 
they do not show any prominent morphological charac- 
teristics aside from a brightening at the western end of 



the remnant. One must be cautious to interpret those re- 
sults, however, as the line emissions of Ne and Mg fall in 
complicated parts of the spectrum, and the fitted abun- 
dances of these elements may be correlated with param- 
eters such as the ionization age, which shows similar dis- 
tribution patterns in the west, and the column density, 
which is very high in that region. It is clear, however, 
that Ne and Mg both show a strikingly different mor- 
phology to Si, S, and Ar, or to Fe, and are much more 
similar to each other than to any of the other elements. 

2.3.2. Forward Shocked Regions 

Broadly speaking, the distinctions between ejecta- and 
FS-dominated regions in Cas A are readily apparent, 
with differences in temperature, ionization age, and el- 
ement abundances. To carry out a survey of the ejecta 
mass, however, we must either model both the ejecta and 
CSM components for each spectrum or else accurately 
identify the specific regions where the reverse-shocked 
ejecta make the dominant contribution to the emission. 
Given the scope of the spectral analysis, we have adopted 
the latter approach. Multi-component fits can be diffi- 
cult to constrain reliably, particularly if one of the com- 
ponents is relatively weak, and thus would require more 
individual attention than is feasible for a sample of thou- 
sands. Consequently, our next aim is to identify and 
eliminate regions whose spectra can be completely asso- 
ciated with the forward shock. 

We evaluate the presence of thermal emission associ- 
ated with the forward shock by fitting a second set of 
plane-parallel shock models to every spectrum, but this 
time with element abundances appropriate for the CSM. 
The optically emitting quasi-stationary flocculi (QSFs) 
in Cas A are understood to be circumstellar mass loss 
from the progenitor. While abundance measurements 
for QSFs are limited to a small number of knots, these 
generally show an order of magnitude enhancement of N 
and sometimes also of He (|Chevalier fe Kirshnerl [19781) . 
Theoretical calculatio ns for t h e pre supernova composi- 
tion are also given by lArnettl ((1996), where the models 
allow the elements H, He, and N, all apparently present 
in Cas A, to exist simultaneously at a narrow tempera- 
ture range near log T (10 9 K) = -1.5 (their Figure 7.6). 
At that temperature the abundance of He is 3 times the 
solar value relative to H by number, and that of N about 
15 times the solar value. As these abundances for He and 
N are broadly consistent with the observational measure- 
ments, we proceed to adopt them for our fits, along with 
solar values for the remaining elements, as representative 
CSM element abundances. 

About 1209 regions gave reasonably good fitting results 
(x 2 < 1-2) with the vpshock model and these QSF ele- 
ment abundances, and are thus assigned to the forward 
shock. They are distributed mainly in the remnant's 
outer rim and southwest interior, as would be expected 
based on the 4-6 keV X-ray contin uum image that high - 
lights the forward-shocked regions dGotthelf et al.l[200l . 
Their average temperature is 2.2 keV, and their ioniza- 
tion ages are rather narrowly distributed with an average 
value of 2 x 10 11 cm~ 3 s. The se values can be ass e ssed i n 
the context of the models of lLaming fe Hwang I {2003), 
which give the current density of the CSM at the forward 
shock at about 1.5-2 cm -3 . Considering the r~ 2 depen- 
dence of the circumstellar density, the forward shock will 
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have encountered much denser material in the past and 
the present-day ionization state of the forward shocked 
material is expected to be relatively advanced. The mod- 
els give values of the ionization age in the 10 11 cm~ 3 s 
range, approaching 10 12 cm _3 s; they also indicate that 
gas is rather hot, with temperatures from 2.5-4 keV. This 
is reasonably close to the average values of the tempera- 
ture and ionization age that we find in our region, though 
the fitted spectra don't show as broad a range in ioniza- 
tion age as is predicted. 

To the forward-shock regions identified solely by the 
thermal emission model above, we must also add those 
that have a strong nonthermal contribution. To identify 
these, we devise a rough diagnostic for the smoothness of 
the X-ray spectrum. We bin each background-subtracted 
spectrum at each significant line feature and continuum 
interval (some of these cover only a narrow energy range) , 
compute the ratio of counts for each major line feature 
relative to counts in an adjacent continuum bin, and take 
the sum of these ratios for all the line features. The dis- 
tribution of this quantity for all the spectra has two over- 
lapping peaks; we take the spectra associated with the 
lower peak (corresponding to weak lines in the spectrum) 
and perform further fits with a composite plane-parallel 
pshock plus power-law model. For a cutoff in x 2 < 1-2 
for these fits, we associate 206 additional regions with the 
forward shock, giving a total of 1415 regions with spectra 
that are consistent with emission that can be associated 
with the forward shock alone. 

Figure 3 summarizes the distribution and spectral 
characteristics of the forward shock regions. Their lo- 
cations echo the 4-6 keV continuum maps s hown by 
iGotthelf et ail (|200ll) and IHwang et all pOOl . As al- 
ready noted, their temperatures and most especially their 
ionization ages are rather narrowly distributed, more so 
than for the sample as a whole. From the Figure, the 
peaks in the distributions correspond roughly to kT = 
2.2 keV and log (n c t) = 11.25, or n c t = 1.8ell cuT 3 s. 
These regions will not be considered further here as we 
focus our study on the ejecta hereafter. We will un- 
dertake a detailed consideration of the thermal emission 
associated with the forward shock in subsequent work. 

2.4. Spectral Survey of the Ejecta and the Presence of 

Pure Fe 

We associate with ejecta the remaining, more than 
4000, regions that are inconsistent with forward shocked 
material. For these regions, the basic spectral model is 
the simple one-component plane-parallel shock that has 
already been presented and discussed. As noted above, 
we take the view that, for the ejecta sample as a whole, 
it is justifiable to assume sufficient ejecta dominance to 
neglect the forward shock component. Representative 
examples of all the various types of spectra seen in Cas 
A are shown in Figure 5. 

For some ejecta regions, the one-component spectral 
model is clearly inadequate to describe the ejecta emis- 
sion in that it fails to accou nt for the Fe K blend. Thi s 
has already been noted by IHwang fe Laming! (2009). 
While there may be significant deficiencies in the atomic 
data that are used for the thermal emission models, par- 
ticularly in the nonequilibrium ionization case, these can- 
not explain so large an effect as we see. In some cases, a 
strong Fe K blend is completely unaccounted for by the 



model, which otherwise characterizes the spectrum well. 

Given the strong chemical inhomogeneities seen and 
expected in core-collapse supernova remnants, it is plau- 
sible that Fe-rich ejecta could be superposed along a 
given line of sight with ejecta of more normal composi- 
tion, and that each could have distinct plasma conditions 
as well as abundances. Moreover, nucleosynthesis models 
predict the formation of nearly pure Fe ejecta by a-rich 
freeze-out during complete Si burning; this is in addition 
to ejecta composed of a range of elements from Si and Fe 
fo rmed by incomplet e Si-bur ning. 

IHwang fe Laming I ()2003l ) and IHwang et all (j2004j ) 
have reported the presence of highly enriched Fe ejecta 
in the southeastern region of Cas A, with one region par- 
ticularly appearing to be a "pure Fe" cloud. We take an- 
other look at this region using the 1 Ms observation. Fig- 
ure 4 shows a portion of the southeastern region of Cas 
A and the roughly 2.5" extraction region that we use for 
this cloud; this is smaller than the region used by Hwang 
et al. 2003. The background is taken locally from the 
vicinity of the source region, and we compared results us- 
ing three different background spectra. The background 
was modeled in this case, rather than subtracted, and the 
model included in the fit to the unbinned, unsubtracted 
source spectrum. The results for all three background 
subtractions gave fitted Fe/Si abundance ratios that are 
enhanced over the solar values by factors of 10-20 by 
number, or factors of 20-40 by mass. The right panel 
of Figure 4 shows the fit giving the highest Fe/Si ratio. 
Note the absence of recognizable emission from S or Ar, 
and the weakness of the Si emission; the emission lines 
at energies above 7 keV are of H-like Fe. These results 
demonstrate that there is emission in Cas A that is plau- 
sibly associated with pure Fe. 

For a subset of 700 regions with the strongest Fe line 
emission and the poorest fits in the Fe K band, we tested 
two options for the added component representing either 
shocked circumstellar medium or pure Fe and Ni ejecta. 
For the former, we take a plane-parallel shock with CSM 
abundances and the average temperature (2.2 keV) and 
ionization age (3 x 10 11 cm _3 s) obtained for the forward 
shock regions, all as described above. For the latter, we 
found that a simple nonequilibrium ionization compo- 
nent, with a single value of the ionization age was suf- 
ficient, with the temperature fixed at 1.95 keV and the 
ionization age fixed at 8 x 10 11 cm~ 3 s. These values 
are representative of the mo st enriched Fe ejecta spec - 
tra, such as those studied by IHwang fe Laming I (2003); 
the fits, however, are not strongly sensitive to the exact 
values used so long as they are of this approximate mag- 
nitude. We fix the parameters of the second component 
whenever possible because the limited spectral resolution 
makes it very difficult, even at these exceptionally high 
levels of signal-to-noise, to constrain both model compo- 
nents independently. As before, errors were calculated 
for key parameters to help ensure finding the true mini- 
mum in x 2 - Only one of these 700 ejecta+CSM fits gave 
X 2 < 1 in the 6-7 keV Fe K region, leading us to conclude 
that the presence of a forward shock component can not 
explain the strong Fe K emission in these spectra. On 
the other hand, a majority of the vpshock+NEl fits did 
give acceptable x 2 values in the 6-7 band, although the 
improvement was sometimes due to a better match to 
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the continuum rather than the line emission. An ex- 
ample is shown in the final two panel rows of Figure 
5, where the same spectrum is shown with the original 
single-component vpshock fit, a two component vpshock 
fit for normal ejecta and shocked CSM, and finally a two 
ejecta-component fit where the added ejecta component 
is taken to be pure Fe (and Ni). 

On a broader scale, we identify 2982 spectra where a 
composite ejecta model providing additional Fe emission 
may be needed, based on a fitted Fe abundance for the 
basic single-component model that is above 0.3 solar. A 
composite model is implemented in cases where the % 2 
per degree of freedom v for the single-component model 
is greater than a user defined threshold. We focus here- 
after only on two component models where the second 
(NEI) ejecta component corresponds to pure Fe. The 
three panels of Figure 6 show examples of the Fe mass 
distribution in the second, "pure" , Fe component for val- 
ues of the threshold for x 2 /v of 1.1, 1.2, and 1.3. In the 
idealized case that fit errors are normally distributed, the 
model linear, and varianc es known such that v is well de- 
fined (jAndrae et al.|[2010h , these values would correspond 
to probabilites of 0.1, 0.01, and 0.001, respectively, that 
the correct model had been fitted. Independently of the 
selected threshold value, the spatial distribution of the 
added Fe distribution is seen to retain its main morpho- 
logical features. 

3. SNR MODELS 

In this section we provide an overview of our mod- 
els for the evolution of Cas A. These are t r eated in a 
manner similar to that in lLaming fc Hwang I ([2003D and 
iHwang fc Lamingl (|2003f) . for a remnant with an ejecta 
density profile with a constant core and a n=10 power- 
law envelope (after Matzner & McKee 1999), expanding 
into a circumstellar wind (density oc r~ 2 ). We discuss 
modifications here. 

We modify the models to accommodate SNR expansion 
in to a stellar wind "bubble" as presented and discussed 
in IHwang fc Lamingl ([2009) . The chief motivation for 
this involves the interpretation of infra-red light echoes 
associat ed with Cas A as repro cessed shock breakout ra- 
diation (jDwek fc Arendtl I2008D . This radiation is con- 
centrated in the UV-EUV spectral region, and would be 
strongly absorbed in photoionizing any intervening neu- 
tral gas, which presumably would be a relic of the red su- 
pergiant presupernova stellar wind. Most of this opacity 
would lie close to the progenitor, and so a small "bub- 
ble" in the circumstellar medium could allow the break- 
out radiation to escape and illuminate the surrounding 
dust. This is especially true if the exploding progenitor 
were compact, and the duration of breakout radiation 
short, as a longer shock breakout would probably pro- 
duce sufficient photoionization to burn its way through 
the surrounding wind. We note that the lib prototype 
SN 1993 J appears to have had dense wind very close 
to the progenitor sin ce it was detected in X-rays a few 
days after explosion (|Fransson et al.lll996l ). While Cas 
A resembles SN 1993J in many other respects, it would 
appear to differ on this point. 

Another update corrects the motion of the reverse 
shoclfl In previous work we have used a model of 

5 Some other minor modifications and corrections to Appendix 



the ejecta where a uniform density core is surrounded 
by a power law envelope. The evolution wh ile the re- 
verse shock is in the outer envelope is given bv lChevalierl 
(1982); this solution is then coupled to the known asymp- 
totic behavior of the Sedov- Taylor solution. In this 
phase, we have previously held the reverse shock speed 
constant at its value at the envelope-core transition, al- 
though it is known that some small acceleration should 
occur during the core propagation. We now imple- 
ment a reverse shock acceleration in this phase of evolu- 
tion following the r esults given for a suite of models by 
iPatnaude fc Fesenl (|2009h . By matching to their results, 
we write for the reverse shock speed and radius: 



v r 



n—3 l E D '■• 



n—3 Ied t co 



43 



43 



(1) 

*,(2) 



where v r = R r /t — dR r /dt. Here s is the power law in- 
dex for the CSM density profile, n that for the ejecta 
envelope, Rb and Vb the forward shock radius and speed, 
respectively, and Ied the "lead factor" (the ratio of 
forward to reverse shock radii, Rb/R r )- The quantity 



x = (40.7AM ej /pR 2 b ) 



l/(3-a) 



with p being the density of 



the CSM at the forward shock in H atoms (or equivalent 
mass) per cm 3 , t core is the time in years when the reverse 
shock hits the ejecta core, and M e j is the ejecta mass in 
solar masses. 

In Table 1 we give various models with the range of 
plausible ejecta masses 2-4 M , all designed to match 
as far as possible the know n dynamics of Ca s A. We 
take a distance of 3.4 kpc (jReed et al.l I1995D and an 
explosion date of 1671.3 ± 0.9 as determined from the 
proper motions of 17 hi gh velocity outer ejecta knots 
(jThorstensen et ahl 120011 ). These knots are identified on 
plates taken as long ago as 1951, providing the longest 
time baseline for such studie s, and are as s umed to be un- 
decelerated. More recently, Fcs en et al.l (|2006f ) have de- 
termined a best explosion date of 1681 ± 19 from a sam- 
ple of 126 relatively undecelerated knots located along 
the NW limb, observed over a time baseline of 9 months 
from March to December 2004. With 1671 as the earliest 
possible explosion date, the remnant age may be up to 
333 years in 2004. 

Further constraints on the model are provided by the 
observed shock velocities and radii. A number of au- 
thors have made measurements of the expansion rate of 
Cas A, both for its bright ejecta ring (IVink et al.l l"l998: 
iKoraleskv et al.lll998t [DeLanev et al.ll2004f) and the for- 
ward shock (|DeLanev fc Rudnick 2003). T he mo st re- 
cent summary is given bv lPatnaude fc Fes cn (2009]), who 
incorporate measurements over the longest time baseline 
(2000-2007). They find that the expansion is slightly 
slower in the N and NW filaments, consistent with the 
known presence of denser CSM on the western limb of 
Cas A. We take the average expansion 0.31" yr _1 from 
other parts of the remnant, which (at a distance of 3.4 
kpc) gives a forward shock velocity of 5050 km s -1 . In- 
cluding all measurements gives a slightly lower average 
expansion of 0.30" yr -1 , and an inferred forward shock 

A in Laming & Hwang (2003) are collected in Appendix A to this 
paper 
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velocity of 4850 km s -1 . 

The forward shock radius is measured from the 
wispy structur e s seen in 4-6 keV continuum emission. 
IGotthelf et all (piffll use data taken in 2000 and de- 
termine an average radius for the northwest sector of 
153 ± 12", which translates to 2.52 ± 0.2 pc, but the 
true shock radius should probably be determined from 
the outermos t of these of the s e stru ctures rather than 
the average. iHelder fc VmFI (|2008l ) perform a decon- 
volution for various azimuthal angle ranges and give a 
slightly larger forward shock radius of 160" (2.637 pc). 
There is also a small variation in the radius around the 
limb. We take a value of 2.6 pc for the forward shock 
radius in 2004 when the Chandra VLP observation was 
taken. 

The reverse shoc k ra dius is harder to q uantify. 
IGotthelf et all (l200l and IHelder fc Vinkl (2008) give a 
range of r adii 1.52 - 1.73 pc, d epending on location in the 
remnant. iMorse et all (2004) measure the reverse shock 
velocity from proper motions derived from two Hubble 
Space Telescope WFPC2 images separated by two years. 
At a reverse shock radius of 106" - 122" (1.75 - 2.01 pc) 
the reverse shock velocity with respect to the expanding 
unshocked ejecta is of order ~ 2000 km s _1 . 

The close match of Cas A's forward shock expansion 
parameter of 0.66 to the Sedov- Taylor value for expan- 
sion into an inverse square density profile has focussed 
interest in low ejecta mass models that allow faster evo- 
lution to the Sedov- Taylor limit. Even so, models gener- 
ally predict too high a shock velocity and/or too small 
a radius for the forward shock compared to observa- 
tions. To i ncrease the forward shoc k radius, it has been 
suggested (|Hwang fc Laming"ll2009t and see above) that 
Cas A exploded into a small bubble in the circumstel- 
lar medium that developed with a transition from a slow 
dense red supergiant wind to a fast tenuous wind shortly 
before explosion. Another possibility is that cosmic ray 
energy losses at the forward shock provide extra decel- 
eration to sufficiently re duce the expansion parameter 
(|Patnaude fc Fesenll2009f ). 

The expansion parameter for a radiative blast wave can 
be different from the Sedov- Taylor limit in the case that 
radiative losses from the shock are significant, or if the 
shock dissipates energy via cosmic ray acceleration. The 
two cases may be trea ted in a similar manner. Following 
iLiang fc Keiltyl (|2000l ) and generalizing to the case of a 
stellar wind preshock density profile p cx r~ s , the radius 
of a radiative blast varies with time as 

rKt l/[4— (3-.)a] j (3) 

where 7 is the adiabatic index of the gas, 

a = {2 - 7 + V /(2-7) 2 + 4( 7l -l)| /4, 

and 

71 - { a/1 + e ( 7 2 - 1) - l} x 2/ (7 - 1) e - 1 

is the adiabatic index modified by the loss of a frac- 
tion e of the shock energy to radiative (or cosmic ray) 
losses. If e = so that 7 = 71, then a = 1/2 for all 
7 and r oc t 2 ^ 5 ~ s K Note that energy going into cosmic 
rays that remain trapped with their energy in the shock 



does not change the dynamics, except through modify- 
ing 7. If all the shock energy is radiated (e = 1), then 
a = l — 7/2. Taking 7 = 5/3, the expansion parameter 
changes from 0.4 to 0.28 for s = 0, and from 0.67 to 0.55 
for s — 2. Consequently, rather small deviations from the 
unmodified expansion parameter are expected for more 
modest cosmic ray energy losses. Cosmic ray energy 
losses of 10% only change the s — 2 expansion parame- 
ter from 0.67 to 0. 6 5, dec reasing to 0.63 for 30% losses. 
iPatnaude fc Fesenl ([20091 ) find variations in the expan - 
sion parameter of this order. iBedogni fc Messina! (|1984f ) 
give full self-similar solutions for several such cases. 

The p lasma density ah e ad of the forward shock is de- 
rived bv lWillingale et al.l (|2003l) from fits to XMM data. 
The baryon mass and filling factor quoted in their Table 
3 give a preshock density of hydrogen atoms or equivalent 
mass of 1.47 amu cm' 3 , assuming a shock compression of 
a factor of 4. Their quoted value for the mass of shocked 
CSM, 8.31 Mq, combined with the dimensions outlined 
above gives a larger preshock density of 1.99 cm~ 3 . These 
values give the time invariant quantity pr 2 = 10 — 13.5 

cm pc 2 where the density p at the forward shock is in H 
atoms (or equivalent mass) cm -3 and the forward shock 
radius 77 is in pc. 

Similar results are obtained from the fi ts of thermal 
bremsstrahlung spectra to Suzaku data (iMaeda et al.1 
120091 ) . interpreted as emission from shocked circumstellar 
medium. They fit the X-ray continuum with a variety of 
composite models. The preshock ion number density de- 
rived from their thermal bremsstrahlung and power law 
fit is 2.1/ (Z 2 ) cm" 3 , where (Z 2 ) = I>^7E 4 ^ 
is an average charge per ion coming from the thermal 
bremsstrahlung emissivity. For gas composed of equal 
number densities of H and He, {Z 2 ) = 1.7 and the 
preshock density is 1.24 cm -3 in H atoms of equivalent 
mass. Taking fit results using thermal bremsstrahlung 
plus SRCu10 or a cut-off power law gives a higher den- 
sity of 1.5 cm~ 3 and pr 2 — 8.4 — 10 cm _3 pc 2 . 

4. DISCUSSION 

4.1. The Elemental Composition of the Cas A Ejecta 

In previous works (|Laming fc Hwangl 120031 : 
IHwang fc Laming I l200l we have discussed the use 
of models of SNR evolution to interpret the fitted 
ionization age of ejecta spectra in terms of the time 
elapsed since reverse shock passage, hence allowing 
the inference of a Lagrangian mass coordinate. This 
is successful in a limited number of carefully selected 
regions of Cas A, but is evidently not the case in general. 
Figure 2 shows the ionization age value strongly peaked 
near log (n e t) ~ 11.3. Converting ionization age to 
mass coordinate, the implication is that all the observed 
ejecta are piled up at a mass coordinate of 0.2 (where 
marks the center of the ejecta and 1.0 the outermost 
extent), instead of being uniformly distributed with 
mass coordinate. It appears that the ionization age is 
not an accurate indicator of ejecta mass coordinate. 
We suggest that this must arise from the interaction 

6 The SRCut model is a widely available model that provides the 
most maximally curved nonthermal spectrum that is plausible, i.e., 
it imposes an exponential cutoff on th e power-law energy spectrum 
of the synchrotron emitting electrons (Reynolds & Keohanc 1999). 
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of the reverse shocked ejecta with secondary shocks 
propagating within the SNR shell, which "refresh" the 
ionization age. These secondary shocks are a direct 
consequence of density inhomogeneities. The forward 
and reverse shocks that encounter them are both trans- 
mitted through the density structure and also reflected 
back into previously shocked ga s. Such a sc enario 
has previously been considered (|Laming| l2001al lbh as 
a possibility for electron acceleration with a view to 
explaini ng the hard X-ray emission of Cas A. More 
recently, llnoue et al.l (|2010f ) consider a similar model for 
ion acceleration at shocks. 

Taking a typical electron density of 200 cm~ 3 
(|Lazendic et al.l 120061 and see below), an ionization age 
of 2 x 10 11 cm~ 3 s implies a shock interaction around 30 
years prior to the observation. Given this "obscuration" 
of the ionization age determined by the reverse shock, we 
estimate the density in each fit region from the emission 
measure assuming 

EM = ^2 n t n e ZfV = UeV^n.Zf = n\V (Z 2 ) (4) 

i i 

where (Z 2 ) = J2i n iZ 2 / J2i n i z i = J2i niZf/n e , and the 
sum i is over all ions in the plasma. This generalizes 
the usual definition for a hydrogen dominated plasma 
to the case of heavy element rich supernova ejecta, 
where the continuum emission is dominated by thermal 
bremsstrahlung emission going as Z 2 . The plasma vol- 
ume is V, and the plasma electron density then follows 
with an appropriate assumption about V. The simplest 
assumption is that V = A 3 / 2 (where A is the area on the 
sky of the spectral extraction region) , and we modify this 
to V = A 3 / 2 x / where / is the filling factor. The fill- 
ing factor is a needed correction for unresolved density 
structure in the ejecta, particularly for SNR evolution 
into a stellar wind (p oc r~ 2 ), as a strong density spike 
is expected in the ejecta density profile. As discussed in 
Appendix B, we characterize the filling factor / by the 
density scale length L p , where / = L p /y/A, and derive 
an approximate value for L p of 2.5" for a distance to Cas 
A of 3.4 kpc. 

The ejecta mass in each imaged region is then esti- 
mated as EM/n e / (Z 2 ) x £\ riimi/n e , from which the 
masses of individual elements are obtained as fractions of 
the total ejecta mass from the fitted element abundances. 
In Table 2 we give the masses of O, Ne, Mg, Si, S, Ar, 
and Fe resulting from our fits, for assumed filling factors 
corresponding to a density spike thickness of 2.5" and 
5" ; the latter arises as one observes the SNR ejecta shell 
in both front and back along the same line of sight. In 
this latter case, the complete X-ray emitting ejecta mass 
inferred is 2.84 M©, and the total ejecta mass includ- 
ing unshocked ejecta is 3.14 M© with reference to Table 
1; this gives a value of pr 2 from Table 1 at the upper 
end of values inferred elsewhere as discussed in section 
3.1. Extrapolation to other cases is simple, since element 
masses are proportional to \[J . Table 2 also gives plane 
of the sky velocities for each element (calculated from 
the position of the center of mass of each element, and 
assuming homologous expansion), and for the SNR as a 
whole, and compares abundance ratios relative to O with 
solar system values. 



We also comment here that Cas A is tho ught to have 
retained some H at the time of explosion (|Fesenl 120011 : 
iFesen fe Becker! [19911 iChevalier fe Oishill2003[ ) ~and this 
may contribute to the thermal bremsstrahlung emission 
here attributed to O. In this case, our ejecta mass es- 
timate would be an underestimate, and if the H (and 
presumably also He) were not uniformly distributed, the 
recoil velocity we attribute to the O ejecta could be in 
error. It is difficult to assess an uncertainty here, but we 
note that the amount of light elements involved is likely 
rather small and unlikely to alter our conclusion of an 
approximately 3 M© ejecta mass. 

In Table 3 we give both the mass of observed Fe ejecta 
associated with incomplete Si-burning, and that associ- 
ated with the "pure" Fe component which can arise either 
from complete Si-burning or a- rich freeze out. We list the 
results for a range of initial \ 2 l v f° r the single vpshock 
fit, as discussed in section 2.4. We also consider the the 
quantity y(1/2) 1^/2 xU ^ 2 ~ 1 ex P — xdx. In the case that 
the fit parameters are normally distributed, and if addi- 
tional conditions are met as discussed at the end of sus- 
bsection 2.4, this quantity would be the probability for 
degrees of freedom v that the single component fit used 
the correct model. The last column column of the table 
gives the fraction f reg of the 2982 ejecta spectra consid- 
ered that were fitted with the composite vpshock + NEI 
model based on the threshold \ 2 l v cutoff for the single 
component model. As can be seen from the table, de- 
manding a higher threshold for x 2 l v (or lower estimated 
probability) to implement the second "pure" Fe fit com- 
ponent reduces Fe Q . At a cutoff value of x 2 l v = 1.2 cor- 
responding to an estimated probability of 0.01, the frac- 
tion f reg begins to be much greater than the estimated 
probability, arguing that we are seeing a real physical 
effect. 

4.2. The Mass and Distribution of Fe 

Our analysis of more than 4300 ejecta regions gives 
a mass of X-ray emitting Fe in the Cas A ejecta of 
~ 0.09— 0.13M©, depending on the assumed filling factor 
and threshold value of reduced x 2 f° r acceptance of the 
"pure" Fe component. There are uncertainties in this es- 
timate, as it is conceivable that the Fe ejecta are clumped 
in a different manner than the other ejecta, possibly as 
in the simulations of Hungerford et al. (2005), and thus 
require a different filling factor than we have used. If, 
for example the filling factor for Fe should be smaller, 
this would yield a smaller Fe mass. The Fe m ass is com- 
parable to that expected (|Eriksen et al.1 12009). but gives 
only the Fe detected in the reverse shocked X-ray emit- 
ting ejecta. According to our models, an extra 0.18 - 0.3 
M© of ejecta may be unshocked, interior to the reverse 
shock, and visible primarily in the infra-red. 

Emission at the center of the remnant that is believed 
to be from unshocked ejecta is seen in infrared Spitzer 
observations, but the emission is of [Si II]. As discussed 
by DeLaney et al (2010), it is coincident with free- free 
absorption seen in the radio by Kassim et al. (1995), and 
appears to correspond to cool (< 1000 K), low density 
gas that has been photoionized. DeLaney et al. con- 
clude, however, that it is unlikely that a significant frac- 
tion of these unshocked ejecta could be Fe. Relatively 
little infra-red Fe emission is observed at all in the rem- 
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nant, a nd the main [Fe II] line at 26 fj, is blended with 
[O IV]. Ilsensee et all (j2010ft demonstrate that [O IV] is 
more plausibly the dominant component of this blend, 
though they do claim a 2a detection of [Fe II] 17.94 /im 
from the center when all spatial bins are summed. Un- 
blended [Fe II ] emission at 17. 94 /im (jEnnis et al.ll2006ft 
and 1.64 /^m (|Rho et ail 12003ft is more convincingly de- 
tected, but appears to be associated primarily with the 
bright ejecta ring, and not the interior. [Fc II] lines are 
also absent or very weak in optic al and near infra-red ob - 
servations of Cas A ejecta knots. iHurford fc Fesenl (1996) 
only report marginal detections of [Fe I I] 8617 A in spec- 
tra of fast moving knots (FMKs), and iGerardv fc Fesenl 
(|2001h report similarly on the [Fe II] lines between 1 and 
2 /im in FMKs, although these transitions are easily de- 
tected in the spectra of quasi-stationary flocculi (QSFs). 
lEriksen et all (2009) detect [Fe II] 1.257 /mi and 1.644 
/jm lines from 15 FMKs and 4 QSFs, and use these to 
estimate the reddening. 

Further consideration of the dust in Cas A sheds light 
on the possible presence of unshocked Fe. Large masses 
of cold dust had be en reported in Cas A from submil- 
limeter observations ([Dunne et al.ll200"3l ). but it has been 
pointed out that some of the detecte d emission originates 
from foreground mole cular clouds (|Krause et al l 1200^ : 
Wil son fc Batrlall2005h . This emission can also be ex- 
plained with < 10 -3 M Q of conducting dust n eedles (i.e., 
Fe) th at are formed in the ejecta (Dwek 2004). Rh o et all 
( 2008) favor FeO dust to explain the Spitzer infrared 
spectra and conclude that Fe dust may be present at 
masses up to lQ ~ 2 Mp,. 

More recently, INozawa et al.l ([20100 suggest the pres- 
ence of significant quantities of cool dust, a conclu- 
sion that is s upported by recent infrared observations 
with Herschel (iBarlow et al.l l2010ft. AKARI and BLAST 
(|Sibthorpe et all l2010D .~The cool dust has a temper- 
ature of about 35 K based on the infrared flux densi- 
ties, and is presumably unshocked ejecta as it is confined 
to the central regions of the remnant. The emission is 
consistent with a silicate dust composition, and the in- 
ferred dust mass is about 0.06-0.075 M Q . Nozawa et al. 
(2010) 's calculations are for dust formation in a Type lib 
event (with an eye toward Cas A) and indicate that little 
dust would be associated with the innermost Fe-Ni layer 
due to the extended radioactive heating: the gas density 
drops too low before the temperatures are low enough 
for Fe or Ni grains to condense. Up to about 1O _3 M 
of Fe could be locked up in FeS grains, but highly Fe- 
rich ejecta would not necessarily have a dust signature, 
as appears to be borne out by the infrared observations. 
Chcrchneff & Dwckj (|2010l ) take a chemical kinetic ap- 
proach and find rather more FeS grain formation in a 
20 Mq model SN with unmixed ejecta (for Cas A, the 
most appropriate of the cases they consider) , with about 
0.021 Mq of Fe in FeS. A mass this high is unlikely for 
Cas A, however, given that its highly stripped progenitor 
underwent a Type lib event. Dust will form relatively 
more efficiently in a normal Type II event because the 
overlying stellar envelope restricts ejecta expansion and 
preserves higher densities favorable for dust condensa- 
tion. 

While Fe appears to be scarce both in the shocked 
ejecta and the unshocked infrared emitting ejecta at the 



center of the SNR, Fe may yet be present in the un- 
shocked ejecta given that the intrinsic emissivities of the 
infra-red lines of Fe II are rather weak compared to, say, 
Si II. Another complication is that there may be a popu- 
lation of low-density Fe ejecta, whose composition and 
distribution are currently unconstrained. Sufficiently 
large clumps of Fe might be expected to reside in low 
density regions due to inflation by radioactivity. Erik- 
sen (2009) note that the lack of detectable Fe lines in the 
Spitzer data indicate that any Fe material still interior to 
the reverse shock must indeed be at a low density, though 
they also suggest that the absence of Fe is more likely to 
be an abundance effect. While an accurate assessment of 
infra-red limits must be deferred to another work, strong 
outward mixing of Fe into at least the He layer has been 
inferred for the Type lib prototype SN 1993J b ased on 
modelling its light curve ([Shigevama et al]|1994D . Cas A 
resembles SN 1993 J, not only in its explosion spectrum 
(seen as a light echo, Krause et al. 2004), but in having a 
similar and low pre-explosion main sequence mass, and a 
binary companion implicated by its strong pre-supernova 
mass loss. 

If our tentative conclusions about the paucity of Fe in 
the unshocked ejecta survive more detailed scrutiny, then 
almost all of the Fe ejected by the supernova is now well 
outside the reverse shock and visible in X-rays, with very 
little left in the center of the remnant. This would be 
true even for the "pure" Fe component, which could be 
the ashes of a-rich freezeout. The inferenc e that the Cas 
A ne utron star has a carbon atmosphere (|Ho fc Heinkd 
2009) would seem to support the notion that the inte- 
rior ejecta composition is dominated by ele ments lighter 
than Fe. Along these lines, recent work bv lOuved et al.l 
(2011) raises the interesting possibility that under the 
right conditions, a quark nova of the incipient neutron 
star could produce 44 Ti by spallation at the expense of 
56 Ni, and also cause C formation consistent with the in- 
ferred composition of the compact object in Cas AQ 

Closely related to explosive Fe nucleosynthesis is the 
production of 44 Ti. Cas A is one of only a handful of su- 
pernova remnants showing clear evidence for 44 Ti in the 
ejecta. An initial mass of 44 Ti of 1.6 x 1O _4 M0 has been 
inferred for Cas A from observations wit h the Comp- 
ton G amma Ray Observatory/COMPTEL (jlvudin et all 
[1994ft . Be ppoSAX and the INTEGRAL IBIS/IS GRI in- 
strument (jVink et al.N2001t IRenaud et al.l 12001 . Based 
on spherically symmetr ic simulations, designed to match 
the C a s A progenitor (fYoung et all 2006; Er iksen et al.l 
12009ft . IMagkotsios eFall pOlOh predict ~ 1O" 4 M of 
44 Ti, and 0.25 M Q of 56 Ni. These are both within a 
factor of 1.6 of measured values (taking the 5" filling fac- 
tor in Table 2), but are discrepant in opposite directions 
so that the mass ratio 56 Ni/ 44 Ti~ 2500 compared with 
the observed ratio of ~ 10 3 . A s imilar outcome of the 
earlier cal culations oflThe et al.l (11998ft and o t hers (re- 
viewed in lThe et all 12006ft led lNagataki et al.l (|1998ft to 
suggest that this ratio may be reduced in an asymmet- 
rical explosion. Another way that 44 Ti production could 
be enh anced relative to F e is a quark nova, as discussed 
above (|Ouved et al.H2011ft . 

7 In this scenario, the compact object is predicted to be a radio- 
quiet quark star surrounded by material rich in C. 
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The results of Magkot sios et al.l (|2010ft also enable an 
estimate of the Fe expected from a-rich freeze out com- 
pared with that formed in incomplete Si-burning. Their 
Figure 17 suggests that approximately half the Fe should 
be associated with a mass ratio 56 Fe/ 28 Si> 10 and half 
with 56 Fe/ 28 Si< 10. If we take a \ 2 /v cutoff of 1.2 from 
our Table 3, we infer that approximately 1/3 of the Fe 
in Cas A could have formed in a-rich freeze out and thus 
be in regions with very small 28 Si masses, in qualita- 
tive agreement with the model estimate. In Cas A, the 
"pure" Fe is placed at least 0.75 Mq out from the cen- 
ter in 3 Mq of ejecta by the the charac teristic ionization 
age of > 2 x 10 11 cm~ 3 s. According to iMagkotsio s et al.l 
( 2010), this is well outside the ejecta region where the Fe 
was formed, and into the outer layers where the compo- 
sition would originally have been dominated by O, Si, C, 
and Ne. The fact that essentially all of the Fe so far is 
found at such locations, and not in the center of the rem- 
nant inside the reverse shock suggests the operation of 
a strong instability, sim ilar to that inferred in SN 1993J 
(jShigevama et al J 11991} ). 

Key information about the explosion of Cas A will be 
provided by the distribution of 44 Ti, a nd this will be 
revea led by observations with NUSTAR (Harrison et al. 
I2010D . If a quark nova did take place during the super- 
nova explosion that produced Cas A, the 44 Ti distribu- 
tion should have a strong central component. Unshocked 
Fe in the center of the remnant will also contribute to 
such a component. In the absence of a quark nova, the 
44 Ti is expected to follow that of the pure Fe, which we 
observe so far in X-rays well outside the reverse shock. 
It is interesting to note that in the case of SN 1987A, the 
ejecta are currently heated by the decay chain of 44 Ti, 
which appears to be located in the center of the rem- 
nant , with ballistic ex pansion velocity of order 3000 km 
s _1 (IKiaer et al.H2010Q. The youn gest Galactic remnant, 
Gl.9+0.3 (jBorkowski et al.ll2010D . is more like Cas A in 
this respect, in that most of its 44 Ti in an X-ray bright 
region beyond the reverse shock, and a smaller portion 
unshocked interior. Considering the young age and fast 
expansion of Gl.9+0.3, this may mean that 44 Ti has been 
mixed out even further into the ejecta than is the case 
in Cas A, although here the SN origin of Gl.9+0.3 is 
unknown. 

Some insight into the motion of the "pure" Fe compo- 
nent can be g ained from the motion of 44 Ti observed wit h 
INTEGRAL (jMartin fe Vinkl [20081: iMartin et al.|[2009h . 
The SPI instrument measures the width and the shift 
of the 1157 keV decay line of 44 Ca, which is the final 
step in the chain of decays from 44 Ti. It is found that 
the 44 Ti is receding from the observer with a velocity of 
about 500 km s , and has an intrinsic width, or velocity 
dispersion, of 430 km s . Since the expansion velocity 
of the Fe is close to 4000 km s -1 , if the Fe is associated 
with the 44 Ti, it must be expanding in a direction close 
to the plane of the sky. The velocity dispersion indi- 
cates that the 44 Ti presumably occupies a volume signif- 
icantly smaller than that outlined by the "pure" Fe, with 
a half angle subtended at the explosion center of order 
sin" 1 (430/4000) ~ 6 degrees. 



4.3. The Neutron Star Kick 



The neutron star is known to be recoiling from the 
explosion center inferred from the optical knot motions 
with a velocity in the pl ane of the sky of 330 km s _1 
(jThorstensen et all I2001D . The direction of recoil is 
also nearly perpendicular from the NE-SW axis estab- 
lished by the jet and counter-jet, which presumably de- 
lineate an axis of rotation in the progenitor. This makes 
a kick arising from the action of a standing accretion 
shock instability appealing. It is known, at least in 
the limit of slow rotation, where distortion of the shock 
front from spherical symmetry can be neglected, that 
the i nstability grows fastest in the prograde spiral mode 
(e.g. iBlondin fe Mezzacappall2007t lYamazaki fe Foglizzol 
l2008t llwakami et al.l 120091) 7 i.e., perpendicular to the 
symmetry axis. The pulsar kick imparted by the sec- 
ond supernova that formed the double pulsar system 
PSR B1913+16 is also inferred to have b een largely 
perpe ndicular to the progenitor spin axis (jWex et al.l 
2000). The magnitude of the observed recoil veloc- 
ity is also well within the accessible range generated 
by numerical simulations of such instabilities, albeit 
these are so far for initially nonrotating supernova c ores 
(jWongwathanarat et aHl201ttlNordhaus et al.ll2010D . At 
higher rotat ion rates, it appears that the instability is 
suppressed (|Burrows et al.l 120071) . though the deforma- 
tion of the accretion shock by rot ation may seed an ax - 
isymmetric quadrupole oscillation (Iwak ami et al.1120091) . 

In the case that the kick involves instability at the 
accretion shock, following th e suspicion advanced by 
IWongwathanarat et al.l (|2010f ). we would expect the Fe 
ejecta to recoil predominantly in the opposite direction 
to the neutron star. In Cas A, this appears not to be the 
case. Most of the Fe, and especially most of the "pure" 
Fe, is recoiling towards the east, within 90 degrees of the 
neutron star motion. However, the remnant as a whole, 
whose mass is dominated by O, is actually recoiling to 
the north, close to 150 degrees from the neutron star mo- 
tion, with a plane of the sky velocity of about 700 km 
s _1 , whic h perhaps suggests a hydrody namic origin for 
the kick. IWongwa thanarat et al.l (|2010f i discuss a 15 Mq 
explosion, whereas the Cas A progenitor was consider- 
ably less massive upon explosion such that reverse shocks 
and associated nucleosynthesis during the explosion are 
less likely to occur. Our estimate of the remnant recoil 
assumes homologous expansion, and the approximations 
made in assessing the ejecta mass have been detailed in 
sections 3 and 4.1. 

5. CONCLUSIONS 

We consider our efforts here to be a first attempt at 
a comprehensive view of the X-ray emitting ejecta in 
Cas A on few arcsecond angular scales. We have car- 
ried out at considerable effort the most detailed data 
analysis that was practicable for us, but there is ample 
room for improvement. The most obvious shortcoming 
of this work is the need for a more sophisticated def- 
inition of spectral extraction regions, and hand-in-hand 
with that, yet more extraction regions. In previous work, 
we had individually chosen regions to correspond to iden- 
tifiable knotty or filamentary features. Here we were 
obliged for the sake of (relative) expediency to instead 
use regions that, while adjusted crudely for overall sur- 
face brightness, have been defined without reference to 
coherent structures in the surface brightness. Full use of 
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the Chandra angular resolution in fact requires a more 
sophisticated definition of the extraction regions. The 
main advantage to be obtained, however, is that this may 
allow the use of simpler NEI models with a single ioniza- 
tion age, rather than parallel shock models with a linear 
distribution of ionization ages, and thus a much clearer 
connection with analytical hydrodynamical models (see 
for ex ample, lLaming fe Hwang | |2003; Hw ang fe Laming! 
2003). It may also eliminate the requirement for the 
filling factor correction used here. Finally, the analysis 
would also benefit from a more sophisticated treatment 
of the spectral background. 

We have arrived at a favored ejecta mass of about 3 
Mq based both on the observed SNR dynamics, and 
on an accounting of the X-ray emission and element 
abundances throughout the remnant. The increased sig- 
nal/noise afforded by the long exposure time, and by the 
relatively large size of extraction regions has allowed us 
to identify regions of the remnant where the Fe ejecta 
is to a large extent "pure", in that the mass ratio Fe/Si 
> 10. We seek to interpret this ejecta component as 
alpha -rich freeze out ashes, where the recently discov- 
ered (jlvudin et al.lll994J : iVink et aLll2001l : iRenaud et al.l 
20 06) 44 Ti and its decay products likely reside. Essen- 
tially all the Fe detected so far, and presumably also the 
44 Ti and its decay products, are not located in the cen- 
ter of the remnant but have been ejected well out into 
the overlying ejecta, probably by the action of hydrody- 
namic instabilities during the explosion. Such outward 
mixing of 56 Ni has been seen in the Type lib prototype 



SN 1993 J (|Shigevama et al.lll994D . With respect to the 
distribution of 44 Ti, SN 1987A presents a stark contrast 
to Cas A, with ejecta containing 44 Ti located at the cen- 
ter with a characterist ic expansion speed of ~ 3000 km 
s _1 (|Kiaer et al.l|2010|). Gl.9+0.3 may represent an in- 
termediate case ()Borkowski et al.|[2010D . 

Finally, Cas A has a well documented displace- 
ment of the neutron s tar formed in the explosion 
(jThorstensen et al~l 120011 ). which is close to perpendic- 
ular to the NE-SW "jet" -axis. While con trary to expec- 
tations (e.g. iWongwathanarat et aLll2010l ). the Fe in Cas 
A does not recoil in the opposite direction to the neu- 
tron star, as far as we can tell, but the remnant as a 
whole does. We take this as reinforcing the notion that 
neutron star kicks in core-collapse supernova explosions 
may derive from purel y hydrodynamic mech anisms, as 
recently suggested (e.g. iNordhaus et aLl l2010). Together 
with motion and location of the Fe ejecta, and the like- 
lihood that what is observed as a NE-SW "jet" really 
is an artefact of an asymmetric explosion, we speculate 
that an instability of the standing accretion shock (SASI) 
in a rotating progenitor is the most promising model for 
understanding the explosion that formed Cas A. 
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Petre, and the referee for helpful comments on the pa- 
per. We gratefully acknowledge support through NASA 
grants to the Chandra Guest Observer Program and 
NNH04ZSS001N to the Long Term Space Astrophysics 
Program. JML was also supported by basic research 
funds of the Office of Naval Research. 



APPENDIX 

A. REVISIONS TO APPENDIX A IN LAMING & HWANG (2003) 

We collect here a few minor corrections and modifications to the text in Appendix A of lLaming fe Hwang I (|2003| ). 
Equations Al and A2 should read 

t = 473.6Af e 5 / 6 £ 51 1/ V 1/ V, (Al) 

x Q = 3A3M^ 3 p- 1/3 pc. (A2) 

The expressions for the lead factor and the pressure ratio between forward and reverse shocks following equation A5 
are now more accurately given by 

Zed = 1-0+ 4 + ^- (A3) 

<fe D = (0.65 - cxp (-n/4)) y/l - s/8. (A4) 

Finally, as pointed o ut to us by Elisabetta Micelotta, some of the discussion following equation A8 in 
lLaming fe Hwang I (120031 ) requires clarification. Equation A5 can be derived from the envelope form of equation A8, 
i.e. from equations A7 and A8 in ITruelove fe McKe"e1 (119991) . This procedure require s v eA = \/ (10/3) (n — 5) / (n — 3) 
going from equation 74 to 75 in ITruelove fe McKee I (|1999l ). In their equation 31, ITruelove fe McKee I (|1999t ) give 
Vcore = v (10/3) (n— 5) / (n — 3), which implies w core = v core /v e j = 1. However their deriva t ion of equation 31 as- 
sumes w core — !• in their derivation. A "cleaner" route to equation 75 of ITruelove fe McKee I (|1999l ) is to take their 
values oi Vrnre. and the eje cta mass fraction in the core derived under the condition w core — s- 0, and substitute into equa- 
tion 3 of iChevalierl (|1982t ). The condition i?t, = l-ET>v core t core still gives the time t core when the re yerse shock hits the 
ejecta, core. The only other change connected with this discussion is that equation All for t conn in lLaming fe Hwang I 
(2003) should depend on v core and not v e j. 

In actual fact, the difference between the approximations w core — and w core = 1 is not all that large; v core — 
y/(10/3) (n — 5) / (n — 3) — > y/2 in these limits, which coincide for n — 8. The core fraction of the ejecta mass varies 
between (n — 3) jn and 1, and only v e j which goes from oo to v core between these two limits is sensitive. 
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B. THE EJECTA FILLING FACTOR 

A correction is needed for unresolved density structure in the ejecta. Models of supernova remnants expanding into 
a stellar wind (p cx r~ 2 ) density profile generally show a strong density "spike" in the radial direction at the contact 
discontinuity, and if the reverse shock is still propagating through the outer powe r law portion of the ejecta density 
profile, the density a t the contact discontinui ty is formally infinite (Chevalier 1982(1 . Later in its evolution, the density 
spike is still present (Chev alier fc Ois hi 2003) , with a width corresponding to about 2 arcsec. This is comparable to or 
smaller than the size of many of our fit regions, so we expect that ejecta structures will not always be resolved, and a 
filling factor correction is therefore needed when calculating the plasma electron density (see subsection 4.1). We take 
the volume of each spectral region o be A 3 ! 2 where A is the area imaged in arcsec 2 . This is modified to A 3 / 2 /, where 
/ is the "filling factor" of the ejecta density profile. 

We estimate a quantitative filling factor / to be applied to our fit results in terms of the density scale length L p as 

/ = Lp/y/AB In deriving L p , we follow IHamilton fc Sarazinl ((1984), who give a similarity solution for ejecta structure 
close to the contact discontinuity. Treating the expansion of uniform density ejecta, they use a coordinate system 
with z = at r — v exp t, the outer limit of the ejecta at time t if they were freely expanding without interacting with 
the reverse shock. The location of the reverse shock is z s , giving dimensionless Eulerian and Lagrangian coordinates 
C = z/z s and £ = Zq/z s respectively, where Zq is the initial value of z in a Lagrangian plasma element at the time that 
the plasma was shocked, i.e. z s at the time of shock passage. With these definitions, close to the contact discontinuity 
they find (their equation 13) 

3/5 -6/5 



/ \ 3 



P = Pl 

C=^fe) ^ + Cc (Bl) 

where v = 23/15 for a stellar wind external density profile, £ c is the value of £ at the contact discontinuity, and 
V s = 1/12 is the dimensionless reverse shock velocity. The pressure is p, and p s is the pressure at the reverse shock. 
Neglecting the dependence of p on r in this region, we write 

dp = dp d( d( = 6 p f p \ 3/5 = 6 p p 

dr 9Co d( dr 5z s V s \pJ ^° hz s v (( - Q L p { ' 

so that the density scale length L p = 23 (£ — Cc) z s /18. We calculate the average of this quantity (L p ) — 
£° pL p r 2 dr/ pr 2 dr using 

"'' pr'dr = £ (*) p sC ^r 2 dr = Ps (£) (3/5 " 2/23) m 18/23 J- r 2 (r c - r^dr 



L pP r 2 dr= 2 -^ 



' Ir r I fe) Ps^r 2 dr = Ps (*) (3/5 " 12/23) m 18/23 i J- r 2 (r c - rf'^d, (B3) 



, 3 / 5 (C _ r \ „ . / „ \ (3/5-12/23) 

"767 

The integrals are evaluated using the substitution u — r c — r, with the final result 

{L , 23 f j (1 - r s ,r c f^ - g (1 - ,„/^) 51/23 + g d - r s /r c )^ \ 



18 I f(!" r s /r c f 23 -1(1- r s /r c ) 28 / 23 + f (1 - r s /r c ) 51 / 23 
Including the variation of p, L p in equation B5 is modified to 

Lp = z s ((- Co) - ^—7^-7- > ■ (B5) 



dlnp/dln(o 

From Table 1 of Ha milton fc Sarazinl (fl98l . 51np/91nCo varies from -1 at the reverse shock to about -0.6 at the contact 
discontinuity, which reduces the value for L p derived assuming <91np/<91n£o = by 5% to 7%. For completeness, if 
the external density cx r~ 3 , then v = 6/5 and L p = z s (C, — Cc) = r c — r, independent of <91np/<91n(o- The filling factor 
is then calculated as / = L p /VA, with L p and A here expressed in arcsec and arcsec 2 respectively. 

As we discuss in the paper, our data suggest that many secondary shocks reflect from the blast wave and reverse 
shock as they encounter density inhomogeneities in Cas A. We therefore use equation B4 with the tacit assumption of 
constant pressure near the contact d iscontinuity. iTruelove fc McKee I (|1999T j point out that this approximation is also 
made bv lHamilton fc Sarazinl (|1984f ) in their energy equation. With r s /r c — 0.9 we find L p = 0.021r c = 1.26 x 10 17 cm, 
which at a distance of 3.4 kpc s ubtends an angle of 2 . 48 arc sec. In this last step we are extrapolating slightly out of the 
formal range of validity of the IHamilton fc Sarazinl {1984) model. There, z c / z s — (v exp t — r c ) / (v exp — r s t) = 0.953, 
which requires an expansion velocity of the outermost ejecta greater than 20,000 km s _1 for values of r c = 1.95 pc 
and r s = 1.7 pc appropriate for Cas A. It is also worth bearing in mind other caveats to the application of this model. 

8 We are considering emission from a "sheet" of plasma, so only < „ ... , , . 

.. . . b 1 j, r 1 i- where all three dimensions would require correction, 

one dimension is compressed, rather than irom a clumpy medium, 
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Hamilton fc Sarazinl (|1984 ) assume uniform density ejecta, while we adopt the prescription of iTruelove &i McKee I 
(1999) of a core-envelope density profile. We have also introduced the circumstellar "bubble" into which Cas A is 
assumed to have exploded. 
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Hwang & Laming 



TABLE 1 
Model Parameters 



M ej /M @ 










Ftp 


prj 


Rbub 




M ej , x /M @ 


m 


age 




2.0 


1.7 


5376 


2.6 


2611 


1.65 


10 


0.26 


4900 


1.84 


0.694 


328.7 


10086 


2.5 


2.0 


5330 


2.6 


2509 


1.67 


12 


0.28 


4924 


2.28 


0.697 


332.6 


9785 


3.0 


2.35 


5357 


2.6 


2465 


1.70 


14 


0.30 


5012 


2.71 


0.699 


331.6 


9683 


3.5 


2.67 


5353 


2.6 


2416 


1.73 


16 


0.33 


5099 


3.13 


0.699 


331.9 


9555 


4.0 


3.0 


5357 


2.6 


2384 


1.75 


18 


0.35 


5163 


3.55 


0.699 


331.8 


9474 



Note. — Table columns give in order: M e j ejecta mass (Mq), E51 explosion energy (10 51 ergs), vi, forward shock velocity (km s 1 ), 
Rl, forward shock radius (pc), v r reverse shock velocity (km s^ 1 ), R r reverse shock radius (pc), pr 2 density times forward shock radius 
at the current position of the forward shock in H atoms (or equivalent mass)x pc 2 , Rb u b circumstellar bubble radius (pc), v exp ejecta 
expansion velocity at reverse shock (km s -1 ), M e j }X mass of X-ray emitting ejecta (Mq), m expansion parameter, SNR age (years), 
and v C ore expansion velo city of ejecta core-envelope boundary (km s — 1 ). The slope of the ejecta envelope is taken to be n=10 from 
IMatzner fc McKeel ltT999Tl . 



TABLE 2 
Element Masses, Velocities 



element 


M ett2 .B"/M @ 


M elt5 „/M Q 


v x 

(km s _1 ) 


Vy 

(km s _1 ) 


Mel /Mq 


M el /M Q \ solar 





1.80 


2.55 


-176 


706 


1.00 


1.00 


Ne 


0.027 


0.038 


1232 


486 


0.015 


0.22 


Mg 


0.0070 


0.0099 


247 


323 


0.0039 


0.12 


Si 


0.038 


0.054 


-672 


677 


0.021 


0.12 


S 


0.020 


0.028 


-685 


787 


0.011 


0.053 


Ar 


0.010 


0.015 


-656 


920 


0.0056 


0.013 


Fe Sl 


0.075 


0.11 


-486 


420 


0.041 


0.23 


Fe Q 


0.023 


0.033 


-1900 


-104 


0.016 




M x 


2.00 


2.84 


-210 


680 






Utot 


2.18 


3.14 


-210 


680 







Note. — Solar abundance ratios from Grevesse et al. (2010). Inferred masses scale as the square root of the filling factor. The Fe 
masses are taken for a \ 2 P er degree of freedom threshold of 1.2 as discussed in the text and given in Table 3. 



TABLE 3 

FE S j, Fe q Masses 



Fe Si,2.5" 


Fe Q , 2 .5" 


Fe Si,5" 


Fe n, 


X> 


1 




freg 


0.061 


0.092 


0.087 


0.13 






10 u 


1.0 


0.072 


0.038 


0.10 


0.054 


1.0 




0.5 


0.29 


0.074 


0.029 


0.10 


0.041 


1.1 




10- 1 


0.21 


0.075 


0.023 


0.11 


0.033 


1.2 




10~ 2 


0.17 


0.075 


0.021 


0.11 


0.030 


1.3 




10 -3 


0.15 


0.076 


0.018 


0.11 


0.026 


1.37 




10- 4 


0.13 


0.076 


0.016 


0.11 


0022 


1.41 




10" 5 


0.11 


0.076 


0.014 


0.11 


0.020 


1.46 




lO" 6 


0.10 


0.076 


0.013 


0.11 


0.018 


1.56 




10- 7 


0.087 


0.076 


0.011 


0.11 


0.016 


1.60 




10- 8 


0.078 



Note. — Fe masses in Mq from incomplete Si-burning (Fegi) and from complete Si-burning and/or ct-rich freeze out (Fe a ), for two 
values of the volume filling factor (corresponding to length scales of 2.5" and 5") and a number of threshold values for the initial single 
vpshock fit. Column 5 gives this threshold value of \ 2 /v, and column 6 the quantity V {v/2) /2 xV ^ 2 ~ 1 ex P —%dx. This last quantity gives 
the probability for degrees of freedom v that the initial fit used the correct model, but only with strict requirements that are not met here 
(see the text). The last column gives the fraction f reg of the 2982 regions considered that were fitted with the composite vpshock + NEI 
model used in the Fe mass calculation. 
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18 Hwang & Laming 




Fig. 1. — Partial set of fitted parameters for the single vpshock model fits. Top row: Maps of \ 2 per degree of freedom (left), and of 
fitted column density Njj in 10 21 cm 2 (right). Second row: Maps of temperature kT in keV and ionization age n c t in cm -3 s (sqrt scale). 
Third row: Maps of fitted Si and Fe element abundances, relative to the solar values of Anders & Grevesse (1989) by number, both on sqrt 
intensity scale, truncated at the high end. Bottom rows: Fitted abundances of the remaining elements, Ne, Mg, S, Ar. The intensities are 
square-root scale and truncated at the high end for S and Ar. 
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Fig. 2. — The distribution of temperature kT and ionization age n e t amongst the 6202 spectral regions fitted with p shock models (see 
Figure 1) : (top) two-dimensional distribution; (lower left) histogram for temperature kT; (lower right) histogram for ionization age n e t 
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Hwang & Laming 





Fig. 3. — (top) Contours showing the 1415 regions that have been associated with the forward shock, and hence excluded for the ejecta 
mass calculation, based on the merged QSF and QSF+PL fits. For these forward shock regions is shown the distribution of fitted (lower 
left) temperature kT (the scale is truncated above 10 keV), and (lower right) logio ionization age n e t. 
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Energy (keV) 

Fig. 4. — (Top) Southeast region of Cas A showing the extraction region used here for the pure Fe cloud of Hwang k. Laming (2003). 
(Bottom) Fitted spectrum, in which the local background has been modelled rather than subtracted and the unbinned spectrum fitted 
using C-statistics. The spectrum is shown binned for clarity, with the black trace giving the total source plus background spectrum, and 
the red trace giving only the source spectrum. 
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Fig. 5. — Sample spectra exemplifying typical spectral types seen in Cas A (first two rows) Spectra associated with the forward shock, 
showing various mixtures of thermal emission associated with shocked circumstellar medium and nonthermal emission, (third row) Typical 
"Si-dominated" and Fe-dominated ejecta spectra, (last two rows) For the same spectrum, a comparison of single-component vpshock ejecta 
models, ejecta plus shocked CSM models, and a model with two ejecta components corresponding to "normal" ejecta and "pure Fe" ejecta. 
For reference, this spectrum had a x 2 / u = 1-98 f° r the single component vpshock model fit. 
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Fig. 6. — The distribution of the "pure" Fe ejecta component in two-component vpshock + NEI model fits where the x 2 /^ cutoff thresholds 
for a single-component vpshock model fit were 1.1, 1.2, and 1.3 (left to right). In the left panel, we show fiducial circles showing the radius 
and position of the forward shock (2.6 pc at 3.4 kpc distance, in yellow) and the contact discontinuity (1.9 pc, in white), and the velocity 
vectors for the NS (white), the "pure" Fe ejecta (light blue), and total ejecta (y ellow) from Table 2 (see subsection 4.2). The circles and 
velocity vectors are centered at the explosion center of Thorstcnscn ct al. (2001). 



